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Abstract 

We study the generalized Forchheimer flows of slightly compressible fluids in heterogeneous 
porous media. The media’s porosity and coefficients of the Forchheimer equation are functions 
of the spatial variables. The partial differential equation for the pressure is degenerate in its 
gradient and can be both singular and degenerate in the spatial variables. Suitable weighted 
Lebesgue norms for the pressure, its gradient and time derivative are estimated. The continuous 
dependence on the initial and boundary data is established for the pressure and its gradient with 
respect to those corresponding norms. Asymptotic estimates are derived even for unbounded 
boundary data as time tends to infinity. 


1 Introduction and preliminaries 

Forchheimer equations [51IB] are commonly used in place of Darcy’s Law to describe the fluid dynam¬ 
ics in porous media when the Reynolds number is large. Their nonlinear structure as opposed to the 
linear Darcy’s equation requires new mathematical investigations. For more thorough introduction 
to Forchheimer flows and their generalizations, the reader is referred to mi, see also [2lll3lll4lll8j . 

In previous articles either for incompressible fluids, e.g. |15[I19] . or compressible ones, e.g. mm- 
the porons media considered are always homogeneous. In reality, however, the porous media 
such as soil, geological media or multi-layer media are not homogeneous. The current paper is to 
start our investigation of Forchheimer fluid flows in heterogeneous porous media. We will develop the 
model and analyze it mathematically. This lays the foundation for onr subseqnent study including 
maximum estimates, higher integrability of the gradient, as well as the structural stability. Such 
analysis will be needed in mathematical theory of homogenization and npscale computation for 
non-Darcy fluid flows in heterogeneous porous media. 

Let a porous medium be modeled as a bounded domain U in space M”' with C^-boundary 
T = dU. Throughout this paper, n > 2 even though for physics problems n = 2 or 3. Let x G ffi” 
and t G M be the spatial and time variables. The porosity of this heterogeneous medium is denoted 
by (/> = </>(x) which depends on the location x and has values in (0,1]. 

For a flnid flow in the media, we denote the velocity by v{x,t) G M”, pressure by p{x,t) G M 
and density by p{x,t) G M+ = [0, oo). 

A generalized Forchheimer equation for heterogeneous porous media is 

g{x,\v\)v =-Vp, (1.1) 


1 


2 


E. Celik and L. Hoang 


where g{x,s) > 0 is a function defined on U x M+. It is a generalization of Darcy and Forchheimer 

equations mm- 

The dependence of function g in (|l.ll) on the spatial variable x is used to model the heterogeneous 
media. For homogeneous media, g is independent of x. For instance, when 

g{x,s) = a, a + Ps, a + f5s + , a + 'yrnS^~^, (1-2) 

where a, /3, 7 , m € (1, 2], 7 ^ are empirical constants, we have Darcy’s law, Forchheimer’s two term, 
three term and power laws, respectively, for homogeneous media, see e.g. mm- Moreover, many 
models of two-term Forchheimer law obtained from experiments, see [2], have a and fd in m 
depending on the porosity (j). For heterogeneous porous media, (p = 4>{x), thus, these coefficients 
become functions of x. This motivates the x-dependent model (|l.ll) . 

In this paper, we study the model when the function g in (jl.ip is a generalized polynomial with 
non-negative coefficients. More precisely, the function g is of the form 

5 (x, s) = ao(x)s"°-|-ai(x)s"^aAr(x)s“^, s > 0, (1.3) 

where N > l,ao = 0 < ai < ••• < are fixed real numbers, the coefficient functions ai(a;), 
a 2 {x), ..., aN-i{x) are non-negative, and ao{x),a]\f{x) > 0. The number on is the degree of g and 
is denoted by deg( 5 ). Such a model ()1.3p is sufficiently general to cover most examples in [2]. 

From (|l.lh one can solve for v in terms of Vp and obtain the equation 

V = —K{x,\Vp\)Vp, (1.4) 

where the function K -. U x —)• is defined by 

K{x,^) = — -^^ for X e D, ^ > 0, (1.5) 

g{x,s{x,()) 

with s = s{x,^) being the unique non-negative solution of sg{x, s) = 

Equation (II.4p can be seen as a nonlinear generalization of Darcy’s equation. The case when 
aj(x)’s are independent of x was studied in depth in dJElEKloHIl]. 

In addition to (jl.lj) we have the equation of continuity 

<P^ + V-{pv) = 0, (1.6) 

and the equation of state which, for (isothermal) slightly compressible fluids, is 

—= w, where the compressibility w = const. > 0 . (1-7) 

pdp 

Substituting (11.41) and (11.61) into (11.71) we obtain a scalar partial differential equation (PDE) for 
the pressure: 

1 

|Vp|)Vp) + K{x,\Vp\)\Vp\\ (1.8) 

at w 

On the right hand side of (jl.Sp the constant w is very small for most slightly compressible fluids 
in porous media, hence we neglect its second term and study the following reduced equation 

0(x)^ = -V-(i^(x,|Vp|)Vp). 

(JL 'UJ 


(1.9) 
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(This simplification is used commonly in petroleum engineering. For a full treatment without such 
a simplihcation, see [3].) 

For our mathematical study of (11.91) below, by scaling the time variable t —)■ we can 

assume, without loss of generality, that w = 1. 

Throughout the paper, function g{x, s) in (II.3p is fixed, hence so is K{x,^). The initial boundary 
value problem (IBVP) of our interest is 

f dv 

4>{x)-^ = V ■ {K{x,\Vp\)Vp) ont/x(0,oo), 

' p = 'll) on F X (0, oo), (1-10) 

p{x,t)) = pq{x) on U, 

where po(3;) are 'ijj{x,t) are given initial and boundary data. 

The main goals of this paper are to estimate the solution of (|1.10p in different norms and 
to establish its continuous dependence on the initial and boundary data. Regarding the PDF of 
(ll.lOp . the fact that (j){x) can be close to zero, alone, makes its left-hand side degenerate. In 
addition, as we will see in Lemma o below and the discussion right after that, the K{x, lVp|) is 
degenerate when |Vp| is large, and can be either very small or very large at different x. Therefore, 
we have to deal with a parabolic equation with different types of degeneracy and singularity. For the 
degeneracy/singularity in x, we use appropriate weighted Lebesgue and Sobolev norms. To identify 
the weight functions for these norms, we carefully examine the structure of the PDF in (II.IOD . 
particularly, the function K(x,^). It turns out that the porosity function (p{x) and the function 
Wi{x), which will be computed explicitly in (I1.17p . are the essential weights. In order to derive 
differential inequalities for the weighted norms, we make use of a suitable two-weight Poincare- 
Sobolev inequality, see (12.2p . We then proceed and obtain the estimates for the solution in section 
[21 for its gradient in section [3l and for its time derivative in section 01 The continuous dependence 
in corresponding weighted norms for both solution and its gradient are obtained in section [5l The 
results for large time are particularly emphasized to show the long time dynamics of the problem. 
Their formulations are made simpler than those in the previous works [8l llOI[T2] . 

For the remainder of this section, we present main properties of K(x,^). First, we recall some 
elementary inequalities that will be needed. Let x,y >0, then 

{x + y)^ < {x^ + yP) for all p > 1, (l-H) 

x^ <1 + x"^ for all 7 > /3 > 0. (1-12) 

The following exponent will be used throughout in our calculations 

a = ^^£(0,1). (1,13) 

We have from Lemmas IIL5 and IIL9 in [1] that 

-aiF(x,0<e^^^<0 V^>0. (1.14) 

This implies K(x,^) is decreasing in hence 


1 _ 1 
g{x,0) ao{x)' 


K{x,C) < K{x,0) 


(1.15) 
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The function K{x,^) can also be estimated from above and below in terms of ^ and coefficient 
functions aj(x)’s as follows. Let us define the main weight functions 


M (x) = max{aj(x) : j = 0,... ,N}, m{x) = min{ao(x), aAr(x)}, 

w r \ (^n{xY , w ^ ^ NM{x) 

= , and W 2 {x) = 


2NM{x 

Lemma 1.1. For ^ > 0, one has 

2Wiix) 


+ aN{xY 


<K{x,0 < 


m{x)aN{xY “ 

W2ix) 




and, consequently, 


aN{x) 


wi{x)e-^ - < K{x,i)e < W2{x)e-^. 


Proof. Let s = s{x,(,) be defined in (ll.Sp . Then 

^ = sg{x, s) = ao{x)s + + • • • + a]sfix)s°‘^~^^ > ai\f{x)s°‘^~^^, 


(1.16) 

(1.17) 

(1.18) 

(1.19) 


hence 

Since the exponents aj are increasing in j, then by ( 11 . 121 ) . one has 
1,..., — 1. Thus, we have 

g{x, s) < M(x)(l + s + ... + s“^) < M{x)N{l + s“^). 

Combining (11.201) and (|1.21l) yields 


( 1 . 20 ) 

< 1 _|_ gO-N fQj. j = 

( 1 . 21 ) 


g{x, s) < M{x)N 


1 + 


aN{x) 


- 

Q^Ar + l 


= NM{x) 


1 + 


aN{x) 


NM{x){Y + aN{xY) 
on^xY 


Therefore 

K(x £) = ^ 

g{x,s) NM{x){f,°-+ ONixY) ^“ + aAr(a;)“’ 

which proves the hrst inequality in (jl.lSp . 

Now, using pi.211) 

= sg{x,s) < s ■ M{x)N{l + (1-22) 

Note that g{x, s) > m{x){l + s“^), then by (ll.22j) and pi.20p . we have 

1 

^ ^ ^ m{x)f, 1 ^ m{x)f, aAr(x)“^+i m(x)aAr(a:)^““^“ 

- NM{x) ' 1 - NM{x) ' ^ NM{x) ' 

Therefore, 

g(x,s) m(x)ajv(x)^~°-^'^ 
hence we obtain the second inequality of (ll.lSp . 
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Next, multiplying (ll.lSp by we have 

2Wiix)e 


+ aN{xY 


<K{x,i)e <W2{x)i 


2—a 


(1.23) 


The second inequality of (jl.23p is exactly that of (ll.lhp . For the hrst inequality of (jl.lhp . if > 
a 7 v(x) then (|1.23l) gives 


Thus, for all .^ > 0 

K{x,OC^ > Wiix){e-‘^ - UNixf-n = Wiix)^-^ - Wi{x)aN{xY-^ 

Note that 

VFi(x)aAr(x)^"“ 


aN{x) 


^ aN{x) ^ a^ix) 


(1.24) 


(1.25) 


(1.26) 


2NM{x) - 2N - 2 

Hence (jl.23p . (I1.25P and (|1.26l) yield the first inequality of (11.191) . The proof is complete. □ 

We now discuss some characters of the PDF in (jl.lOp . On the left-hand side, the porosity 
(f){x) can be close to zero, hence giving the degeneracy in variable x. On the right-hand side, the 
dependence of K{x,^) on ^ as seen in (jl.lSp shows that the PDF is degenerate in |Vp| as |Vp| —)> oo. 
Moreover, since the weights Wi{x) and W 2 {x) can tend to either zero or infinity at different location 
X, then, thanks to (I1.18h again, so can K{x,^). Therefore the PDF can become singular and/or 
degenerate in x. The above fact about the weights Wi{x) and W 2 {x) is supported by practical 
models in [2]. For example, the two-term Forchheimer law (i.e. iV = 1) has the coefficients oq and 
ai going to zero as (/>—)• 1, and to infinity as ^ 0. In this case, (/> is required to be in (0,1). For 

heterogeneous media, constant (f) becomes function (j) = Yix) £ (0,1) and it may be close 1 or 0 at 
different values of x. Therefore, thanks to the mentioned behavior of ao(x) and ai(x), the weights 
Wi{x) and W 2 {x) possess the stated property. 

Lastly for this section, we recall an important monotonicity property for the PDF in (ll.lOj) . 
Lemma 1.2 (c.f. [T], Proposition III.6 and Lemma IIL9). For any y,y' G M”', one has 
iK{x,\y\)y - K{x, \y'\)y') ■ {y' -y)>{l- a)K(x, max{|?/|, \y'\})\y - y'^. 

In order to estimate the pressure gradient, similar to mm. we will make use of the function 

re 

lo 

Same as (96) of [T], we have the comparison 

K{x, Of < H{x, 0 < 2K{x, Of- (1.28) 

Combining (ll.28p with (11.191) gives 


re 

H{x,0 = / K{x,\fs)ds for X G [/, 0. 

Jo 


(1.27) 


VFi(x)e"-“ - 


aAr(x) 


<H{x,0<‘^W2ix)Y 


2—a 


(1.29) 
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2 Estimates for the pressure 

We start analyzing the IB VP (ll.lOp . To deal with the non-homogeneous boundary condition, let 
^{x,t) be an extension (in x) of 'ijj{x,t) from boundary T to U. Our results are stated in terms of 
T, but can be easily converted to i/j, see e.g. [7]. 

Let p = p — ^, then we have 

dp 

(l){x)^ = V ■{K{x,\Vp\)Vp)-(t){x)^t on[/x(0,oo), (2.1) 

p = 0 on r X (0, oo). 


The analysis will make use of the following two-weight Poincare-Sobolev inequality 


'u 


\up4>{x)dx ) < cp { I Wi{x)\Vuf °'dx 


1 

2 —a 


( 2 . 2 ) 


for functions u in certain classes that satisfy u = 0 on P. 

For some classes of functions (/>, IFi, u such that the inequality (12.211 is valid, see e.g. mm- 
Here we give a simple example that (12.2p holds under the so-called Strict Degree Condition 


deg(fi() < 


n — 2 


(SDC) 


Note that in the three dimensional cases (n=3), (SDC) reads deg( 5 () < 3, hence it holds for the 
commonly used two-term, three-term and power Forchheimer models, see m- 

We recall the standard Sobolev-Poincare’s inequality. Let IF^’'^(17) be the space of functions in 
with vanishing trace on the boundary. If 1 < < n then 

< c||V/||M(t/) for ah / e (2.3) 


where the constant c depends on q, n and the domain U, and q* = nq/{n — q). 

One can easily verify that (SDC) is equivalent to 2 < (2 — a)*. 

Assume (SDC). Let q < 2 — a such that r = q* > 2. Let c be the positive constant in (12.3p . 
Assume further that 


def 

Cp = c 


u 


Wi{x) ^-°--‘idx 


2 — a — q 
(2-a)q 


(f){x)^-^dx 


u 


r-2 

2r 


< OO. 


(2.4) 


Let u € Then by Holder’s inequality and standard Sobolev-Poincare inequality (12.3p 


^ J \u\'^(j){x)dx 



(j){x) ’’-2 dx 


u 


< c 


\Vu\'^dx 


4>{x) ^-^dx 


u 


r-2 

2r 


Since q < 2 — a, applying Holder’s inequality again to the second to last integral, we obtain (12. 2 p 
with Cp defined by ( 12 .4p . 

The above example shows the validity of (12.2p for reasonable cp, Wi while u belongs to a standard 
Sobolev space. Nonetheless, for the purpose of this paper, it suffices to assume, without focusing 
on technical weighted Sobolev spaces, that the inequality (|2.2p always holds true for u = p, as well 
as u = P in section [5l 

Notation. The following notations will be used throughout the paper. 
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If /(x) > 0 is a function on U, then define 

LP{U) = |■u(x) : \\u\\lp^u) (^ J^f{x)\u{x)\Pdx^ 


\ i/p 


< oo 


(2.5) 


Notation || • \\]^p will be used as a short form of || • WiP^uy 

• We will use the symbol C to denote a generic positive constant which may change its values 
from place by place, may depend on the domain U, dimension n and the Sobolev constant cp 
in ( 12 . 21 ) . but not on individual functions aj(x)’s, and not on the initial data and boundary data. 
Constants Cq,Ci,C 2 , ■ ■ ■ have fixed values within a proof, while di,d 2 , ■ ■ ■ are fixed positive constants 
throughout the paper. 

• The notation pt stands for Similarly, pi^t = P 2 ,t = etc. 

• For a function f{x,t), we denote f{t) = 

In this section, we derive estimates for p{x,t) in 
Lemma 2.1. Ift>0 then 

^ J p^{x,t)4'{x)dx + j iF(x, lVp(x,t)|)|Vp(x,t)p(ix < C'Go(t), ( 2 . 6 ) 

where C > 0 and 


Go(t) = + J ao{x) ^|VT(x,t)p(ix + J Wi{x)\V^{x,t)\‘^ ‘^dx + (^ J \Ft{x,t)\‘^4>{x)dx 


with 


= 


/ a]\f{x)dx. 

Ju 


2-a 

2{l-a) 

(2.7) 

( 2 . 8 ) 


Proof. Multiplying equation (12.ip by p{x, t), integrating over U, and using the integration by parts, 
we obtain 

/ p^cfdx = — / K{x,\Vp\)Vp -Vpdx — / pFt4>dx. 

2 dt Jp Jp Jp 

Substituting p = p — T into the first integral on the right-hand side, and applying Cauchy-Schwarz 
inequality give 


f p‘^(pdx+ [ K{x,\Vp\)\Vp\'^dx = I K{x,\Vp\)'V'F ■ Vpdx — f p'Ftcfdx 
2 dt Jp Jp Jp Jp 


< 


[ iF(x,|Vp|)|VT||Vp|dx+ [ \p\\^t\Hx = h+h- (2.9) 

Ju Ju 


• For Ii in (12.9p . applying Canchy’s ineqnality, we have 

h= / iF(x,|Vp|)|VT| • |Vp|dx < ^ [ K{x,\Vp\)\Vp\‘^dx+ [ K{x,\Vp\)\V^\^dx. 

Ju 4 Ju Ju 

Estimating the last integral by (ll.lSp . we get 

1 


h < 


[ K{x,\Vp\)\Vp\‘^dx 

Ju 


+ 


'U 


ao{x] 


iVTFdx. 


( 2 . 10 ) 
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• For I 2 in (12.91) . applying Holder’s inequality and using the weighted Sobolev-Poincare inequal¬ 
ity (j2.2p . and applying Young’s inequality with powers 2 — a and we have 


h< 




\P\ 


'U 
< 2 ““^ 




u 


<cp[ I iyi(x)|Vp| °'dx 

u 


1 

2 —a 




u 


lFi(x)|Vp|2-“dx + c( [ 

^ Ju 


Id-* 


2—a 
2{l-a) 


JU ^Ju 

By triangle inequality, (11.111) . and relation (11.191) we estimate 

[ lFi(x)|Vpp-“dx < 2^-“ [ (lFi(x)|Vpp-“ + iyi(x)|V^p-“)dx 

Ju Ju 

(iF(|Vp|)|Vp|2 + 


< 2 ' 


< 2 


1—a 


>U 


iF(|Vp|)|Vppdx + 2-“Hi+ 2 ^-“ / lFi(x)|V^p-“dx. (2.11) 

Ju 


Therefore, 

K{x,\Vp\)\Vp\‘^dx + CBi + C J lFi(x)|VT|2-“dx + c(y' (2.12) 

Combining (12.91) . (12.101) . and (12.121) . we have 


IjL 

2 dt 


[ p‘^(j)dx + ^ [ K{x,\Vp\)\Vp\‘^dx 

lu ^ Ju 


<C [ ^-iVTpdx + CHi+ C [ lTi(a:)|VT|2-“dx + C 

Ju Mx) Ju 



2-a 

2{l-a) 


Thus we obtain (12.Op . 

For the sake of future estimates’ simplicity, we replace Bi in (12.7p by 

= maxjili, 1}. 


□ 


(2.13) 


Thus, (12.61) gives 


^ [ p‘^{x,t)4>{x)dx + I K{x,\'Vp{x,t)\)\'^p{x,t)\'^dx < CG{t), 

dt Ju Ju 


(2.14) 


where 


G{t) = G[mt) = + [ ao(x)-^|VT(x,t)pdx+ [ lFi(x)|VT(x,t)|2-“dx 

Ju Ju 

+ (y \'^t{x,t)\^(l>{x)dx^^^^^. (2.15) 

With this change, we have 


G{t) >1 Vt > 0. 


(2.16) 
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Using (12.1111 and (|2.14l) . we derive 


A 

dt 


'u 


p + 2“ / lUi(a:)|Vp| “dx 


'f/ 


< 


4 ; [ p^^dxF [ K{x,\Vp\)\Vp\'^dx + CB^ + C [ lUi(x)| 
dt Ju Ju Ju 


VM'^-^dx 


<CG{t) + CB^ + C / lUi(x)|V^r"“dx < UG(t), 

Ju 


which, by setting di = 2“ proves 


4 [ p'^{x,t)4>{x)dx + di [ Wi{x)\Vp{x,t)\‘^ °‘dx<CG{t). (2-17) 

dt Ju Ju 


Applying inequality (12.21) to tt = p and utilizing it in ()2.17p give 
d 


dt 


b 


p‘^{x,t)(p{x)dx < —d 2 { I p^{x,t)4>{x)dx] ‘‘ +CG{t), 


u 


2 — 0 , 
2 


(2.18) 


where d 2 = diCp“^. 

This nonlinear differential inequality enables us to obtain estimates for p in terms of initial and 
boundary data. They are described by the following function and numbers. 

Let Ad(t) = be a continuous function on [0,oo) that satisfies 


Denote 


A4(t) is increasing and JA{t) > G{t) Vt > 0. 


A = limsupG(t) and B = limsup[G'(t)] . 

t^OO t^OO 


Note from (j2.16p that 


JA{t) >1 Vt > 0, and ^ > 1. 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


Assumptions. Throughout the paper, we assume that each solution p{x, t) and its corresponding 
function 'L(x,t) have enough regularity in spatial and time variables such that all calculations 
are carried out legitimately. Also, the time dependent quantities such as the above G{t) and others 
introduced in subsequent sections are required to belong to C([0, oo)), and when needed, C^((0, oo)). 
The purpose of this requirement is to allow application of Gronwall’s inequality and other types of 
estimates in Appendix [Al 

Theorem 2.2. 


(i) If t > 0 then 


/ p^{x,t)(t>{x)dx < / p^(x, 0 )(^(x)(ix + C'Ad(t) 2 -a. 

Ju Ju 

(2.22) 

(ii) If A < oo then 


/ 2 ^ 
limsup / p {x,t)'f{x)dx<CA‘^-°-. 

t^oo Ju 

(2.23) 

hi) If B < oo then there 

is T > 0 such that for all t > T 



f p^{x,t)4>{x)dx < C{B^^ + G{t)'^). 

Ju 

(2.24) 










10 


E. Celik and L. Hoang 


Proof, (i) Define y{t) = p'^{x,t)f>ix)dx. We rewrite (I2.18P as 

y'{t) < + CG{t), (2.25) 

2 — ^ 

where = Cqz^^ with Cq = ^2 ■ Using nonlinear Gronwall’s inequality in Lemma lA.lf ii. 

we have for all t > 0 

y{t) < 2/(0) + (p{CM{t)), 

hence obtaining ( 12 . 221 ) . 

(ii) Applying Lemma I A. If hi to the differential inequality (I2.25|) . we get 

2 2 

limsup 2 /(f) < C'limsupG(t) 2 =^ = CA^^, 

t^oo t—>-oo 

which proves (|2.23p . 

(iii) Finally, note that ip{z) < g>o{z) Cq{z + z^) for z > 0, where 1 < 7 = < 2. 

Clearly, (/Jq ^(y) < Lp~^{y), then we have from (I2.25p that 

y'{t) < -ip-^{y{t))pCG{t). 

Hence by Lemma lA.21 there is T > 0 such that for all t > T 

y(t) <G (1 + H^+G(t)^), 

which, together with (I2.16p . yields (I2.24p . □ 


3 Estimates for the pressure’s gradient 

In this section we estimate the weighted L^^“-norm for the gradient of p. Due to the structure of 
equation (j2.ip . we start with estimates for H{x, |Vp(x,t)|) defined by (ll.27jl . and will use relation 
(I1.29P to derive the ones desired. We define 

Gi{t) = Gi[^]it)= [ ao{x)-^\V^t{x,t)\^dx. (3.1) 

Ju 


Theorem 3.1. 

(i) For t > 0, 


/ H{x,\Vp{x,t)\)dx < e G / H{x,\S/p{x,G)\)dx 

Ju Ju 

+ g(^J p^{x,0)4){x)dx + (t) + j e~G^~GGi{T)dT^ . (3.2) 


(ii) If A < 00 then 


limsup / H{x,\Vp{x,t)\)dx < g(a^-°- 

^00 Ju ^ 


+ lim supGi(t) 


t—^OO 


(3.3) 
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Proof, (i) Multiplying equation (I2.ip by pt, integrating over U, and using integration by parts we 
have 


lu 


pjcfdx = — / K{x.,\Vp\)Vp-Vptdx — / pt^t<pdx 

Ju Ju 

/ K{x,\Vp\)'\/p-Vptdx + / K{x,\Vp\)'Fp ■'S/'Ftdx — / pt'^tf'dx 

Ju Ju Ju 

H{x,\Vp\)dx+ / K{x,\Vp\)Vp ■ V'Ftdx — / pt^t4>dx. 

Ju Ju 


'u 
1 d 


‘2^ dt Ju JU 

Let e > 0. Applying Cauchy’s inequality, we derive 

1 d 




1 


Pt4>dx + I H{x,\Vp\)dx < e I K{x,\Vp\)\Vp\'"‘dx + — I K{x,\Vp\)\VFt\'"‘dx 
^ dt Ju Ju 4e Ju 


1 


+ ^ / Pt4>dx + 


1 


\"Ft\4>dx. 


lu 


'u 


By using (jl.28l) to estimate the first term on the right-hand side, and using (ll.lSp to estimate the 
second term on the right-hand side, we obtain 


^ JuPj(pdx + ~ JuH{x,\Vp\)dx 


< e 


lu 


H{x,\Vp\)dx + ^ [ aoix) ^\V'Ft\‘^dx + ^ [ I'Ft\‘^(f)dx. 
Ju ^ Ju 


Thus, 


[ pi(j)dx + ^ [ H(x,\Vp\)dx < 2e [ H(x,\Vp\)dx + —Gi(t) + G(t). 
Ju dt Ju Ju 

From (12.141) and (ll.28p . we have 




'u 


H{x, \Vp\)dx < CGit). 


(3.4) 


(3.5) 


Combining ^3.41) and (13.5p with e = 1/8, we obtain 

■ 7 - [ p^(t>dx + ^ [ HixJVplxjPDdx + [ picl)dx + - [ HlxJVplxjPDdx 

dt Ju dt Ju Ju 4 Ju 

<C{G{t) + Gi{t)). (3.6) 

We rewrite the first term on the left-hand side and apply Cauchy’s inequality as follows 

-/ 

dt Ju 

hence 

H{x,\Vp{x,t)\)dx — f pfcjydx + - f H{x,\Vp{x,t)\)dx < 2 f p^cfdx + C{G{t) + Gi{t)). 

Ju ^ Ju Ju 

(3.7) 

Particularly, neglecting the second integral of the left-hand side reduces (13.71) to 


g2 


p'^(pdx = ‘2 1 pptcfdx > —- f Pt(l)dx — 2 / p^cjdx, 

Ju 2 Ju Ju 


-I 

dt Ju 


f H{x,\Vp{x,t)\)dx + - f H{x,\Vp{x,t)\)dx < 2 f p'^Jidx + C{G{t) + Gift)). (3.8) 
Ju 4 Ju Ju 


jL 

dt JU 
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Using (12.221) to estimate the integral term on the right-hand side, and then properties (12.191) . 
(j2.21[) . we have 

^ |Vp(x,t)|)(ix + \Vp{x,t)\)dx 

<2 [ p‘^{x,0)^dx + CM{t)^+C{G{t) + Gi{t)) <C [ p'^ix,0)^dx + GM{t)^ + GGi{t). 
Ju Ju 

Consequently, by Gronwall’s inequality. 


'u 


H{x,\'Vp{x,t)\)dx < e 4* / H{x,\'Vp{x,0)\)dx 

Ju 

-|- C J e“4h“U^y p‘^(^x,0)4>dx + + Gi{T)^dT. (3.9) 


Since A4(t) < A4(t} for all r € [0,t], estimate ()3.2p follows (|3.9p . 

(ii) Applying Lemma lA.lIf ni to differential inequality (13.81) , and using limit estimate (I2.23P , we 
have 


lim sup 

t—>-oo 



iL(x, |Vp(a;,t)|)dx < Climsup / -|-C'limsup(G(t)-|-Gi(t)) 

t—^oo Ju t^oo 

< G{A''^ -|- .A + limsupGi(t)). 

t—>-oo 


Since A. > 1, by (I2.2ip . we obtain (13.3p . □ 

For large time, we improve the estimates in Theorem 13.II bv establishing inequalities of uniform 
Gronwall-type [T71I20] . 

Lemma 3.2. For t>l, 


[ H{x,\'^p{x,t)\)dx + 1- f f pf{x,T)(l){x)dxdT 
Ju ^ Jt-^ Ju 

<g(^J p'^{x,t — ^)4>ix)dx + {G{t) + Gi{T))dT^. (3.10) 

Proof. The proof follows [8] by using basic differential inequalities (13.4j) and (13.5h . 

Integrating (13.5p from t — 1 to t yields 

f p‘^{x,t)'f>dx + I- f f H{x,\Vp{x,T)\)dxdT < f p'^{x,t — ^)(pdx + G f G{T)dT. 

Ju 2 Jjj Ju Ji-_i 

Neglecting first term on the left-hand side, we have 

f f H{x,\Vp{x,T)\)dxdT < 2 f p'^{x,t — ^)4’dx + G f G{T)dT. (3-11) 

Jt-i Ju Ju Jt-i 

Using (13.4p with e = 1/2 gives 

f p‘f{x,t)4>dx + — f H{x,\Vp{x,t)\)dx < I H{x,\Vp{x,t)\)dx + G{G{t) + Gi{t)). 

Ju dt Ju Ju 
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Let s £ — l,t]. Integrating the previous inequality in time from s to t, we obtain 

n pf{x,T)(j)dxdT + f H{x,\Vp{x,t)\)dx 

j Ju 

< f H{x,\Vp{x,s)\)dx + f f H{x,\Vp{x,T)\)dxdT + C f (G(t) + Gi(r))dr 
Ju Js Ju Js 

< f H{x,\Vp{x,s)\)dx + f f H{x,\Vp{x,T)\)dxdT + C f (G(r) + Gi(r))(i7 

Ju Jt-i Ju Jt-i 


Integrating the last inequality in s from t — 1 to t results in 
rt rt 

Jt-l 


n p‘f{x,T)(l)dxdTds+ f H{x,\Vp{x,t)\)dx 
r Ju 


<2 / f H{x,\Vp{x,T)\)dxdT + C f {G{t) + Gi{T))dT. 

J t—1 Ju J t—1 


Using (|3.11l) to estimate the Hrst term on the right-hand side, we have 

rt rt 

Jt- 

< 


n pf{x,T)(j)dxdTds+ f H{x,\Vp{x,t)\)dx 

J Ju 

4 f p'^{x,t — l-)(t>dx + C f {G{t) + Gi{T))dT. (3.12) 
Ju Jt-i 


rt-l/2 Hr 1 /■* /■ 

p^{x,T)4>dxdTds > / / / p^{x,T)(j)dxdTds = - / / pf{x,T)(l)dxdT. 

Jt-i Jt-l;Ju ^Jt-irJu 


For the first integral on the left-hand side, we observe that 

rt rt r /■i-1/2 pt 

p'f{x,T)4>dxdTds > I 

/ 1—1 Js Ju 

Utilizing this estimate in (I3.12j) . we obtain inequality (j3.10p . 

Combining Lemma 13.21 with Theorem 12.21 results in the following specihc estimates. 
Theorem 3.3. 

(i) Ift>l then 


□ 


'U 


H{x,\Vp{x,t)\)dx < c(^ J p‘^{x,0)4>{x)dx F M-{t)‘2-<i'+ Gi(r)(ir^. (3.13) 


(ii) If A < oo then 

limsup / H{x,\Vp{x,t)\)dx < c(a^^ + lira sup / Gi{T)dT 

t^oo Ju ^ t^oo Jt —1 

(iii) If B < oo then there is T > 1 such that for all t > T, 

J H{x,\Vp{x,t)\)dx <c(b^ +G{t)^ + Gi(r)dr). 


(3.14) 


(3.15) 
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Proof, (i) Combining (13.1011 with estimate (I2.22p and property (|2.19p yields 

f H{x,\'\/p{x,t)\)dx < c( f p'^{x,0)(j){x)dx + M-{t — 1)^^ + f {G{t) + Gi{T))di 

Ju ^Ju Jt-i 

< g(^J P^{x, 0)4>{x)dx + M(t)^ + M(t) + Gi(r)(ir^. 

Then using the fact A4(t) > 1 from (j2.2ip . we obtain (j3.13p . 

(ii) Taking limit superior of (|3.10ll . and using limit estimate (|2.23ll . we have 


limsup / i7(x, |Vp|)(x,t)(ix < C'limsupG(t)2-a + Climsup / [G{t) + Gi{T)]dT. 

t^OO Ju t—j-oo t—>-oo Jt—1 


Note that 


Then 


limsup / G(t) dr < limsupG(t). 
i—>-oo Jt —1 t—>-00 


(3.16) 


t—>CXD Ju 


limsup / H{x,\Vp\){x,t)dx < G[A^-°- + A+ \im.sup / Gi(r)(ir). 


t^oo Jt-l 


Estimate (I3.14P then follows since ^ > 1. 

(hi) Using (I2.24p to estimate the term J^p^{x,t — l)(/>dx term in (I3.inp . we obtain 

J H{x,\Vp{x,t)\)dx < g(^B~ + G{t - 1 )^ + (G(r) + Gi(r))dry (3-17) 


Note from Lemma lA.41 that 


G{t) < G{t) + B + 1 VTG[t-l,t]. 


Hence (j3.17l) implies 

I H{x,\Vp{x,t)\)dx <g(b^+ {G{t)+B + l)^+ {G{t)+B + l) + Gi(T)dT). 


Then inequality (I3.15P follows by using (|l.lip . (|1.12p and the fact G{t) > 1. □ 

Remark 3.4. (a) Compared to (j3.2p . estimate (j3.13ll does not require Vp(x,0). Also, (I3.14p im¬ 
proves (13.31) slightly, particularly when Gi {t) fluctuates strongly in time, (b) The estimate (|3.15p is 
simpler than (I3.17p which is the form usually presented in previous papers PlIUHE]. 

The statements in Theorems o and 13.31 can be rewritten to give estimates for the integral 
JjjWi{x)\Vp{x,t)\'^~'^dx, that is, \\Vp{t)fAj,. 

Corollary 3.5. For t > 0, 


f Wi{x)\Vp{x,t)\‘^ °'dx<e 4* f H{x,\Vp{x,0)\)dx 

Ju Ju 

+ g(^J p‘^{x,0)4>{x)dx + (t) + J e~i'^^~'^'^Gi{T)dT^. (3.18) 
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For t > 1, 

f Wi{x)\Vp{x,t)\‘^~°'dx < c( f p‘^{x,0)(l){x)dx + + f Gi(r)di 

Ju ^ Ju Jt-i 


If A < oo then 


If B < oo then there is T > 1 such that for all t > T, 


lu 


Wi{x)\Vp{x,t)\'^~‘"dx <C{B^ +G{t)^ + I Gi{T)dT 

t-i 


Proof. Using property (ll.29|) . definitions (12.81) and (I2.13P we have 
Wi{x)\Vp{x,t)\‘^~°'dx < I 

Jv 


lu 


(3.19) 


limsup / Wi[x)\Vp{xA)\^ “dx < _)_ iiixLsupGi(t)y (3.20) 

t^oo Ju ^ t^oo ' 


(3.21) 


QaK^) .g H(^x,\Vp{x,t)\) dx < ^ + J^H{x,\Vp{x,t)\)dx. 


Also, from definition ()2.15p . G{t) > With these relations, the estimates (I3.18p . p3.19p . p3.20p . 
(I3.2ip immediately follow (13.21) . (|3.13l) . (13.1411 . (I3.15p . respectively. □ 


4 Estimates for the pressure’s time derivative 

In this section, we estimate the pressure’s time derivative. Let 

q{x,t) = Ptix,t) and q{x,t) = Ptix,t) = Pt{x,t) - 

Then q solves 

da 

0(a^)-^ = V • {K{x, |Vp|)Vp)t - (p{x)^tt on U X (0, oo), 
5 = 0 on r X (0,oo). 

In the following estimates, we use 

def 


G 2 (t) = G 2 [T](t) = / \'Ftt{xA)\ (f{x)dx. 

Ju 


(4.1) 


(4.2) 


Lemma 4.1. One has for t > 0 and e > 0 that 

-f- [ q‘^{x,t)4>{x)dx <—{1 — a) [ K{x,\Vp{x,t)\)\'^q{x,t)\‘^dx 

dt Ju Ju 

+ e [ \q{x,t)\'^(pix)dx + CGi{t) + Ce~^G 2 {t). (4.3) 

Ju 

Proof. Multiplying (|4.ip by q, integrating over U, and using integration by parts we have 

J ^qJ'dx = J V ■ {K{x, \Vp\)Vp)tqdx - j ^ittjjqdx 

= - {K{x, |Vp|)Vp)t • Vqdx - / Fttq4>dx. 

Ju Ju 
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Taking the derivative in t for the first integral on the right-hand side, we derive 




dK{x^\Vp\) (Vp^ Vq) _ r ^(^jVp|)V(7-Vgdx- [ ^t^dx. 

c'? \yp\ Ju Ju 


Using the fact that q{x, t) = q{x, t) — for the first two integrals on the right-hand side, we rewrite 


l_d 

2 dt 


q^4>dx = — 


dK{x,\Vp\)\Vp-Vq\^^^^ f dK{x,\Vp\)Vp-Vq ^^ 


'u 


lu 


di \Vp\ Ju \Vp\ 

— / K{x,\Vp\)\Vq\‘^dx + / K{x,\Vp\)Vq ■ V'^tdx — / ^ttq4>dx. 

Ju Ju Ju 

Next, by derivative property (I1.14p and Cauchy-Schwarz’s inequality we obtain 

[ q‘^J'dx< a [ K{x,\Vp\)\Vq\‘^dx +a [ K{x,\Vp\)\VqWV'^t\dx 
2 at Ju Ju Ju 

- [ K{x,\Vp\)\Vq\^dx + J^K{x,\Vp\)\Vq\\V'^t\dx - [ ^ttq<pdx 


lu 


'U 


<-(l-a) [ K{x,\Vp\)\Vq\‘^dx + {a + l) [ K{x,\Vp\)\Vq\\V'kt\dx + [ \^tt\\q\J)dx. 
Ju Ju Ju 

Let e' > 0. Applying Cauchy’s inequality to the last two integrals gives 

Ju ~ + 1) - (1 - a)) \Vp\)\Vq\^dx 

+ ~rr f K{x,\Vp\)\V'^t\‘^dx + e [ \q\'^4>dx + ^ [ \'ktt\‘^4>dx. 

4e Ju Ju 4 e Ju 

We estimate K{x, |Vp|) in the second integral on the right-hand side by (ll.lSp . then it follows 

1 d 


„ , q'^J’dx < {e'{a + l) - {1 - a)) K{x,\Vp\)\Vq\‘^dx 

2 at Ju Ju 


+ 


tt -j-1 


/ ao{x) ^\Vt\‘^dx + £ / \q\‘^(f)dx- 

Ju Ju 


G2it) 


4e 


Selecting e' = (1 — a)/(2(l -|- a)) gives 


+ CGi{t) + £ j^\q\'^4>dx 


+ 


CG2it) 


which proves ^4.31) . 


□ 


The next theorem contains different estimates of ||pf(t)||L 2 for both small and large time in 
terms of the initial and boundary data. Note that we cannot estimate the norm at t = 0. 

Theorem 4.2. 

(i) For to S (0,1] and t > to, 


[ H{x,\Vp{x,t)\)dx+ f pl{x,t)4){x)dx < gUq^ [ H{x,\Vp{x,0)\) + p‘^{x,0)4>{x) 

lu Ju ^ Ju ^ 


u 

+ J Gi(r)(ir-k A4(t)2^-k J e“3h-^)(Gi(T)-k G 2 (r))dT^. (4.4) 


dx 
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(ii) If t > 1 then 

J pf{x,t)4>{x)dx < c(^J p^{x,0)4>{x)dx + M{t)^ + {Gi{t) + G 2 {T))dT^ . (4.5) 

(iii) If A < oo then 

limsup f pf{x,t)4>{x)dx < c(a^^ + lira sup f {Gi{t) + G 2 {T))dT 

t^co Ju ^ t^oo Jt—1 


t—>oo Ju 

Consequently, 


lims 

t^oo 


(4.6) 


(4.7) 


(4.8) 


isup / pl{x,t)(f{x)dx < C+\inisup{Gi{t) + G 2 

-^00 Ju t^oo 

(iv) If B < 00 then there is T > 0 such that for all t > T, 

f pf{x,t)4>{x)dx < c(b~ + G{t)^ + f {Gi{t) + G 2 {T))d' 

Ju ^ Jt-i 

Proof. Denote I{t) = fjj H{x,\Vp{x,t)\)dx + Jjj q‘^{x,t)(j){x)dx for t > 0. 

(i) Adding ()3.7p to (14.3p with e = 1/4 yields 

7-1 p 

-I{t) + -I{t) < C / p^J)dx + CGsit), where G 3 = G + Gi + G 2 . 
dt 4 

Integrating (13.61) in time from 0 to t, we have 

f I{T)dT < Gi/(0) + Cl [\g{t) + Gi(T))dr 
Jo Jo 

for some Gi > 0. Applying ()4.10l) to t = to) then there exists t* G (0,to) such that 

I{u) < ^ f\G{T)+Gl{T))dT. 

to to Jo 

For t > to, applying Gronwall’s inequality to (14.91) on the interval [G,t], and then combining 
with (14.lip , we have 


(4.9) 


(4.10) 


(4.11) 


I{f)<e 4^* **)/(G) + G [ e 4!* '^1 ( p^{x,T)(l){x)dx + G‘j,{T) 

Ju ^Ju 


dr 


< 


^1(0) + ^ [ (G(r) + Gi(r))dr + G [ e 4!* / p^(x,r)(/(x)dx + G3(r) 

to tQ Jq Jq IJjj 


dr. 


Using ()2.22p to estimate the integral jjjP^{x,T)(f{x)dx yields 

G /■*« 


m<^m + ^ r(G(r) + Gi(r))d7 
to to Jo 


+ C f e 4 !* f p^(x, 0)(/>(x)(ix + A4(r) 2 -a-p ( 73 ( 7 -) 

Jo ^Ju 


dr. 
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Since Ai{t) is increasing, see (|2.19l) . it follows that 

+ ^ [ Gi{T)dT + CM{to) + C [ p‘^{x,0)(j){x)da 
to to Jo Ju 


2 H 1 

+ M{t)^+C e-4(‘-^)G3(r)dr. (4.12) 

Jo 


Also, from (12.211) . JA{t) > 1, then JA{to) < A4(t) < A4(t)2-“. Thus, we obtain estimate (j4.4l) from 

dm. 

(ii) Using (14.3p with e = 1/2, and dropping first integral on the right-hand side, we have 


2 dt 




q^Jidx <]: [ q^Jidx + C2{Gi{t) + G2{t)) 


(4.13) 


lu 


for some G 2 > 0. For s G [t — t], integrating (14.131) in time from s to t gives 

i f q^{x,t)(t>dx f q^{x,s) 4 >dx + ^ f f \q{x,T)\^(l)dxdT + G2 j (Gi(r)-b G2(r))(ir 
j ^ J J s J cy J s 

]- f q'^{x,s)4>dx + 1- f [ \q{x,T)\'^4>dxdT + C 2 f (Gi(t)G 2(r))d7 
^ Ju ^ Jt-\ Ju Jt -1 


< 


Next, integrating in s from f — ^ to f, we have 


1 1 
2 ' 2 


f q'^{x,t)(t>dx < I- f [ q'^{x,s)(l)dxds + ]- - I- f [ \q{x,T)\'^(j)dxdT 

Ju ^ Jt-^ Ju 4! 2 J^_l J^ 

^ [ {Gi{t) + G 2 {T))dxdT = ^ f f q‘^{x,s)(l)dxds + ^ f (Gi(r)-b G 2 (r))dr. 
Jt—l ^ Jt—h Ju ^ Jt—l 


Hence, 


lu 


q'^{x,t)4’dx < 3 f f q'^{x, s)(j)dxds + C f {Gi{t) + G 2 {T))dT. 

Jt-\ Ju Jt-i 

Using ()3.10l) to bound the first integral on the right-hand side, we obtain 

f pf{x,t)4>dx<G f p'^{x,t — ^)4>dx + G f {G{t) + Gi{t) + G 2 {T))dT. (4:.14) 

Ju Ju Jt-i 

Combining ()2.22p and (I4.14p gives 


/ Ptix,t)H- 
Ju 


x<G p'‘(x,0)4>{x)dx + GJA{t — 1)^-°-+ G / (G(t)-|-Gi(t)-|-G 2 (r))dr. 

Ju Jt-i 

Again, by properties (12.191) and (I2.2ip . estimate (14.511 follows. 

(iii) Taking limit superior of (I4.14j) and using (I2.23h . we obtain 

limsup f pf {x,t)4>dx < GA^^ + GA + Glim sup f (Gi(r)-|-G 2 (r))dr, 
t^oo Ju t^oo Jt—l 

which yields (14.61) . The estimate (14.7p follows (14.61) and property (I3.16P for functions Gi and G 2 in 
place of G. 
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(iv) For sufficiently large t, estimating the first term on the right-hand side of (|4.14l) by (12.241) . 
and then applying Lemma I A. 41 to bound G{t — 1) and G(r) in terms of B and G{t), we obtain 


/ pUx,tWx<C{B^+G{t- 1 )^) + C (G(r) + Gi(r) + G 2 (r))dr 

Ju Jt-i 

< CB^ +C{l + B + G{t))^ +C{l + B + G{t)) +C [ (Gi(r) + G 2 (r))dr. 

Jt-i 

Then (|4.8I) follows by simple manipulations using inequalities (11.111) . (11.121) . □ 


5 Continuous dependence 

In this section, we establish the continuous dependence of the solution on the initial and boundary 
data. 

Let pi{x,t) and p2{x,t) be two solutions of (11.91) with boundary data ’ipi{x,t) and 'ip 2 {x,t)^ 
respectively. For i = 1 , 2 , let t) be an extension of ipiix, t), and define pi = Pi — 'hj. Denote 

P = Pi — p2, <1> = 'Ll — ^2 and P = pi — p2 = P — 


Then 

dP 

= V • {K{x, |Vpi|)Vpi - K{x, |Vp2|)Vp2) - on D X (0,oo), 

P = 0 on F X (0, oo). 

The weighted norms of P and $ are related by the following differential inequalities. 
Lemma 5 . 1 . For all t > 0 , one has 

d 


dt 


'u 


P^{x,t)<l){x)dx <—d'ihi{t) 2-a I I Wi{x)\'VP{x,t)\‘^ °‘dx] + CD{t)h2{t)'2, 


<u 


2 

2 —a 


— [ P'^{x,t)(j){x)dx <—dihi{t) 2-a f P‘^[x,t) 4 >{x)dx + CD{t)h2{t)^, 


lu 


'u 


dt 

where ^3,(^4 > 0, 

D{t) = J aQ{x)~^\V^{x,t)\‘^dx + J ao(x)“^| V<L(x, ^^ y \Ft{x,t)\‘^cj){x)dx 

2 r 

hi{t) = Bi + '^ H{x,\Vpi{x,t)\)dx, 

7^1 Ju 

2 /. 

h2{t) = l + '^ / H{x,\Vpi{x,t)\)+p‘f{x,t)(p{x) dx. 


i=i JU - 


(5.1) 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 


Di{t) = [ \Ft{x,t)\‘^ 4 >{x)dx, £>2(1)= [ ao{x) ^\V^{x,t)\‘^dx, 

Ju Ju 

2 2 

hsit) = '^\\piit)\\l2, /i 4 (t) = ^ / H{x,\Vpi{x,t)\)dx. 


2 = 1 


2 = 1 


Proof. We define 
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Multiplying equation (15.111 by P and integrating over U give 

/ P ■ Pt(j)dx = / {V ■ {K{x,\'\/pi\)'\/pi — K{x,\Vp 2 \)Vp 2 ))Pdx — / ^tP(pdx. 

Ju Ju Ju 

Using integration by parts for the first integral on the right-hand side, we have 
~ [ P^(^dx = - [ {K{x,\Vpi\)Vpi-K{x,\Vp 2 \)Vp 2 )-VPdx- [ 

dt Ju Ju Ju 

= - [ {K{x, |Vpi|)Vpi - K{x, |Vp 2 |)Vp 2 ) • (Vpi - Vp 2 )dx 

Ju 

+ / {K{x,\Vpi\)Vpi — K{x,\Vp 2 \)Vp 2 ) ■ V^dx — / ^tPJ^dx. 

Ju Ju 

Applying Lemma 11.21 to the third to last integrand, we obtain 


1 ^ r 
2 dt Ju 


P'^4>dx <—{1 — a) / Ar(x, |Vpi| V |Vp 2 |)|Vpi — Vp 2 pda: 


lu 


+ / (iL(x,|Vpi|)|Vpi|+iL(x,|Vp2|)|Vp2|)|V$|dx+ / m\P\cpdx. 


Above, we use the notation |Vpi| V |Vp 2 | = max{|Vpil, |Vp 2 |}- 
For the first integral on the right-hand side, we note that 


|Vpi - Vp2p = |VP + V^>p > - |V^>|^ 


hence. 


^ d f -2 ,, I - a 

P (bdx < - 

2dt ' ^ - 


'u 


'u 


K{x,\Vpi\W\Vp 2 \)\VP\‘^dx + C [ A:(x,|Vpi| V |Vp 2 |)|V$pdx 

Ju 


+ / {K{x,\Vpi\)\Vpi\ + K{x,\Vp 2 \)\Vp 2 \)\V<^\dx + / \<^t\\P\(Jdx 
Ju Ju 

= + Ch + h + h. 

Consider Ii. Let IC{x,t) = K{x, |Vpi(a;,t)| V \Vp 2 {x,t)\). Then by Holder’s inequality, 

Wi{x)a 


Wi{x)\S/Pg °‘dx < ( I IC{x,t)\VPfdxj Ji , where Ji = 


2 —a 


-dx. 


'U JC{x, t) a 

Applying (|1.18p to bound IC{x,t) from below, and then using (jl.lip . we estimate Ji as 

Ju ^ 21 Fi(x) / 


(5.7) 


(5.8) 


<C{ I Wi{x)aNixy-'^dx + J Wi{x){\Vpi\\/ \Vp 2 \y-‘^dx 


< C Wi{x)aN{xy-'^dx + iyi(x)(|Vpi|^-“ + \Vp2rndxj . 

Then by (|1.26l) and (ll.29p . we have 

Ji <c(^J aN{x)dx + J [H{x, |Vpi)| + H{x, \Vp 2 \)]dx'^ = Chi{t). 
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This and (j5.8l) yield 


Ii = j JC{x,t)\^P\ dx>C Wi{x)\VP\^-^dx] 


2 

2 —a 


(5.9) 


For 12, by using (11.151) 


I 2 <C [ ao{x)-^\V^\‘^dx = CD 2 {t). 
Ju 

For I 3 , applying using Holder’s inequality gives 


(5.10) 


'3< E 

i=l ,2 


^ K{x, iVpiDiVpipdx) ^ K{x, |VKl)|Vd>|2dx) ^}. 


Using ()1.28p for the first integral and (11.151) the second integral, and then applying Cauchy-Schwarz 
inequality, we have 


/3 <a/2 V / H{x,\Vpi\)d: 


i=l ,2 • 


lx ] [ I ao(x) 


Thus, 


I 3 < Ch 4 (t) 21 ) 2 ( 1 ) . 


For I 4 , applying the Holder’s inequality gives 


(5.11) 


I 4 < ||P|L2||$i||^2 < (\\pi\\l 2 + \\p2\\L2)\\^t\\Ll = h3{t)2Dl(t)2. 

(p <p q> <p <p 


(5.12) 


Then combining ()5.7p . ^5.91) . (I5.10p . (15.lip and ^5.121) yields 


2 

2 —a 


^ J P'^pdx < -ds ( / VFi(x)|VP|^-“dx ) 


'u 

+ CD2(t) + Ch4{t)^ D2(t)^ + h3{t)^ Di(t)F 


Hence 


2 

2-a 


— J P‘^{x,t)(l>{x)dx < —d 3 hi(t) 2 -a Wi{x)\VP{x,t)\‘^ “dx ] + 00 ^( 1 ), (5.13) 


where 


We estimate 


D^it) = D2(t) + h4(t)2D2(t)2 + h3{t)2Di{t)2. 


Db(t) < C(D2{t) + D2{t)2 Di{t)2)(l h4(t)2 /i3(t)2) < C D{t)h2{t)2. 


(5.14) 


Therefore, (15.2p follows (|5.13l) and (|5.14l) . Finally, using Poincare-Sobolev’s inequality (12.21) for 
u = P, (|5.2p implies (15.3p . □ 
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To describe more specific estimates, we introduce 

2 2 


Po = / p‘i{x,0)4>{x)dx, ^0 = y] / H{x,\Vpi{x,0)\)dx, 

7^1 du ~1 Ju 


i=l 


and, referring to (I2.15p . (I2.19p . (13.ip . (14.21) . define for t > 0 




2 = 1 
2 


2=1 

2 


Gi(t) = G2{t) = ^G2[T,](t). 

2=1 2=1 

In the following, we show that the L^-norm of P{t) for t > 0 can be bounded by the initial 

difference ||P(0)||j;^2 and the difference between the boundary data expressed by H(t). It means 

4> 

that the solution of (jl.lOl) depends continuously on the initial and boundary data. 

Theorem 5.2. For t > 0, 

\\P{t)\\l2 +C (5.15) 

0 _/q 


where 


A4i(t) = ^0 +'Po + + sup G'i(r). 

Tg[0,t] 


In particular, for any T > 0, 

sup \\Pmh<\\m\\h+CMi{T)l f D{t)dt. 

te[o,r] 0 0 Jq 

Proof. Define y{t) = jjj P‘^{x,t)(j){x)dx. We rewrite (15.3p as 

y'{t) < -dihi{t)~^y{t) + GoD{t)h 2 {t)^. 

By Gronwall’s inequality 

Jo 


(5.16) 


(5.17) 


(5.18) 


Let t > 0. Using definition (15.51) of hi{t) and by applying (13.2p to each p = Pi for z = 1, 2, we 
have 

hi{t) <B^ + e-3‘^0 + CVo + GM^{t) + G f e-^^*-^^Gi{T)dT. 

Jo 


Thus, 


hi{t) < Gil-Lo+ Po + + sup Gi(r)) = CA4i(t). 

rG[0,t] 


(5.19) 


Similarly, by (15.61) . estimates (13.21) and (12.221) . we have 

h 2 {t) < 1 + e-^^o + Po + GM{t)^ +G [ e-^^^-^'>Gi{T)dT, 

Jo 
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which implies 

h2it) < CMi{t). (5.20) 

Therefore, we obtain (j5.15p from (j5.18ji . ()5.19p and ()5.2np . 

Now, let T > 0. Neglecting the exponentials in (j5.15p and noting that M-i{t) < M.i{T) for all 
t € [0,T], we obtain (j5.16l) . □ 

Next, we estimate P{t) when f —)• oo. The estimate is independent of the initial data and only 

depends on the asymptotic behavior of 'I'i(x,f), 'I' 2 (x,f), and as t —)• oo. 

hi{t) = oo, then from (|5.17l) and Lemma lA.ll iil we have 

limsup ||.P(f )||^2 <^limsup—= Climsup(5-21) 

t—^-OO ^4 t—>-00 ^ 2 —a t^OO 

where 

R{t) = h2{t)^hi{t)^^ D{t). (5.22) 

To estimate the last limit in (I5.2ip . we define, referring to (12.201) . (13.11) . and (14.2p . the following 
numbers 


4 = = ^limsupG[Tj](f), ^^limsup[(G[Tj](f))'] , 


2 = 1 


2=1 

2 


£—>-cx) 


2=1 


2=1 


£—>-cx) 


Gi = ^limsupGi[4'j](f), 02 = ^limsupG2[4'i](f). 


£—>•00 


£—>■00 


The asymptotic behavior of <h(x,f) as f —)• oo will be characterized by 


T> = limsupZl(f). 

£^oo 


Denote also that 


Ko 


a 2 + a 

2-a 2 “ 2(2-a)' 


Theorem 5.3. If A and 0i are finite, then 

limsup||P(f)||2 2 < C{A^ +giT°V. 

t^OO ^ 


Proof. Note from (12.231) and ()3.3I) that 

~ 2 

limsup/ii(f),limsup/ 12 (f) <C{A^-°- F Oi) < 00 . 

t—^oo t—^oo 


(5.23) 


(5.24) 


Then hi{t) and /i 2 (f) are bounded on [0, 00 ). Thus, hi{t) ^-°-dt = 00 and, consequently, estimate 
(I5.2ip holds. By (I5.22p and (|5.24l) . 

limsupi?(f) < C{A^^ + ^ 1 ) 2 ( 42 ^ + ^ 1 ) 2^22 = C{A'^ + GiY^V. (5.25) 


£—>•00 


Therefore, (I5.23jl follows this and (I5.2ip . 


□ 
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Now, we focus on the case when the boundary data is unbounded as t —)• oo. 
• If t > 0 then, by (12.221) . 


2 

I 2-a 


/ Pi{x,t)4>{x)dx < C\Vo + M{ty 

i=i 

If t > 1 then, by (13.1,41) . 

^ r 2 

J H{x,\Vpi{x,t)\)dx < c(Po +M{t)^ + G'i(r)dr^, 

and, by (14.5p . 

Y] / plt{x,t)4>ix)dx <c(Po + M{t)^ + / (Gi(r) + G' 2 (T))di 

^ Ju ^ Jt-i 


(5.26) 


(5.27) 


(5.28) 


• In case B < oo, then ,S[^i] and ,S['I' 2 ] are finite. Using estimates (I2.24p . (I3.15p . (14.8p for pi 
with i = 1,2, there is Tq > 0 such that for t > Tq, one has 


2 „ ^ 

/ p‘iix,t)4>{x)dx < c(b^^ + G{t) 


2 

2 —a 


(5.29) 


yy / H{x,\Vpi{x,t)\)dx < + / (5i(r)dr^, (5.30) 

^^ JU Jt —1 

2 /• /-t 

J^plt{x,t)4>{x)dx < g{b~ +G{t)^ + ^{Gi{t) + G 2 {T))dTy (5.31) 

• Assume A. = oo. Then limj_j.oo A4(t) = oo, and for sufficiently large t, one has M(t) > Vo- 
From (j5.27l) . (15.281) . and (15.301) . (I5.3ip . we have for large t that 


V / H{x,\Vpi{x,t)\)dx, V / plPx,t)(j){x)dx < GV{t) 

JU 


(5.32) 


2=1 


where 


U(t) = 


(5.33) 


A4(t)2-“ + JYYi{x) + G 2 {T))dT in general, 

+ G{t)^^ + J^_i{Gi{t) + G 2 {T))dT when B < oo. 

With the above preparations, we are ready to estimate ||.P(t)|| r 2 as t —oo in the case A = oo 

4> 

Theorem 5.4. Assume A = oo. If f^V(t) ^-°-dt = oo, then 


limsup ||P(t )||^2 < Climsup (v{tP°D{t)) 

t—^OO ^ V J 


t^OO 


(5.34) 
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Proof. By (15.291) and (I5.30p . or (|5.26l) and (|5.27l) . we have for large t 

h2{t),hi{t) <CV{t). (5.35) 


Combining this with (I5.2ip . we have 

limsup ||.P(t )||^2 < C limsup{h 2 {t) 2 D{t)) < C'limsup(C(t) 2 i/(t) 2 ^L)(t)) 

t^oo ^ t^oo >-oo 

= Climsup(i/(t)^°L>(t)). 

t—)-oo 


This proves (I5.34jl . □ 

The estimate (|5.34l) can be interpreted as follows. As t —>■ oo, even though V{t) —)• oo, if the 

boundary data’s difference characterized by D{t) decays very fast, it can diminish the growth of 

V{t) and result in ||P(t)||j ;^2 going to zero. 

4> 

Now, we turn to the continuous dependence for the pressure gradient. What the results obtained 
below mean for ||VP(t)|L 2 -a are the same as Theorems 15.21 and 15.41 for ||P(t)||j;^ 2 . 

^Wi 4> 

Theorem 5.5. Let Iq G (0,1]. For t > to > 0, 


\\VP{t)\\%.a < CM2(tr°(e-'^^fo^dr) ^dr||p(Q)||2 + D(tf 

V W 


where 


Moreover, 


M 2 {t) = tQ^iUo +Vo) + + sup (G'i(r) + G' 2 (r)). 

TG[0,t] 


(5.36) 


lim sup ||VP(t)|| 

t—^OO 


T ^ — a _ 


< c 


{A^ + gi + G2f^°v 


+ C{A^ +GiT°V. 


(5.37) 


_1 a 

Proof. Multiplying (15.2p by df , we have 

2 

(y j P‘^{x,t)(p{x)dx + Chi{t)^^D{t)h2{t)2 

<Chi{t)^ [ \P\\Pt\f)dx + CR{t) 

Ju 

< C*/ll(t)2-“ ||Pj||j^2 ||P||^2 + CR{t). 

Applying triangle inequality to ||A||r2 gives 

</> 

||VP(t)||2,_, < C'hi(t)^(||pi,t||^2 + \\p2,thimLl + CR{t). (5.38) 

^Wi 4> <t> <P 

By (lO^ . (l5Tm and (IMP . 

R{t) < CMi{t)lMi{t)^D{t) = CMi{tY°D{t). (5.39) 
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By (I13D, 


5 l,t|lL 2 

4> 


(U) + l|P2,t||L2([/) < C(tQ^ f Gi{T)dT + Iq ^{Ho + Vo) 

J 0 

rt 1 

+ 7W(t)2^+ / e-3(^-^)[G2(r) + Gi(r)](iT)^ <CA^2(t)^. 
Combining estimates (|5.19p . (|5.4Up . (|5.15p . (|5.39p with (I5.38P yields 


(5.40) 


+ 




Estimating the first and last A4i(t) terms on the right-hand side by A4i(t) < M. 2 {t), we obtain 

(1061) . 

-_ 2 _ 

Let Vo = (^ 2 -a Taking limit superior of p5.38p . and using the limit estimates p5.24p . 

P5.23I) . (I5.25P and (|4.7p . we have 

limsup ||V.P(t)||^ 2 -a < C{A'^ + Gi)"^+Gi + G 2 )^V^ + CVo 

t-yoo Wi 

< c{A^ + 01 + g2)"-^vh + c{A^ + gi)^°v, 


hence obtaining (15.371) . 

Finally, we derive the gradient estimates for the case ^ = oo. 
Theorem 5.6. Assume .4. = oo. Let V{t) be defined by p5.33p . Suppose 

roo 

/ V~^^{t)dt = oo and lim (172^^ (t))'= 0. 

Then 


□ 


limsup IIVP(t)||22-a < Climsup V{t)'^^°D{t) +(71imsup V{t)^°D{t) 
t—>-00 t^OQ L -I t—>-00 L 

Proof. By (I5.38p . (I5.35p . (I5.32p and (I5.22p . we have for large t that 

||VP(t)f 2 _<, < CV{t)^°\\P\\L2+CV{t)^°Dit). 

Taking limit superior of the previous inequality yields 


(5.41) 


(5.42) 


limsup IIVP(t)||': 2 -a < Climsup 

t^OO t^OO 


Ju 


P\^4>dx] -|-C'limsup(17(t)''°L)(t)). (5.43) 


t^OO 


Consider first limit on the right-hand side of the p5.43p . Let y(t) = Jjj\P{x,t)\'^(j){x)dx. By 
P5.17p and p5.35p . we have for large t that 

y'{t) < -CiV{t)~^y{t) + C 2 D{t)V{t) + 

We apply Lemmato h{t) = CiV~'^(t), f{t) = C 2 D{t)V{t )2 ^ and g{t) = noticing 

that condition pA.2p is met thanks to (I5.4ip . It follows that 

limsup (V{t)^^°y{t)) < Climsup{Vit)^'^°V{t)^Vit)^D{t)) = Climsup{Vit)^^°D{t)). 


Then inequality (I5.42p follows this and (15.431) . 


□ 
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A Appendix 

We collect here some useful lemmas on solutions of differential inequalities. 

Lemma A.l (c.f. [9], Lemma A.l). Let (j) be a continuous, strictly increasing function from [0, oo) 
onto [0, oo). Suppose y{t) > 0 is a continuous function on [0, oo) such that 

v'{t) < -h{t)(t)~^{y{t)) + f{t) Vt > 0, 


where h{t) > 0 , f{t) > 0 are continuous functions on [ 0 ,oo). 

(i) If JG[[t) is an increasing, continuous function on [0,oo) that satisfies M{t) > f{t)/h(t) for 
all t >0, then 

y{t) < y{0) + 4>{M{t)) Vt > 0. 

(ii) If h{T)dT = oo then 

limsupy(t) < (p 

t—)-C50 

Here, we use the notation Pioo) = oo. 

Lemma A.2 (c.f. [8], Proposition 3.7). Let (p{z) = c(z + z^} for all z > 0, where c > 0 and 
1 < 7 < 2 . Suppose y(t) > 0 is a continuous function on [0, oo) such that 

y'{t) < -(p-^{y{t)) + f{t) yt > 0, 


1- fit) 

limsup—^ 
t^oo a[t) 


where f{t) >t) is a function in (^([O, oo)) n C^((0, oo)). 

Assume /3 limsupj^oo[/'(t)]“ is finite. Then there is T > 0 such that 

y{t) < C{1 + fity) for all t > T, (A.l) 


where C = 3[32(1 + c)] 2 - 7 . 

Proof. We track and calculate the constant C explicitly. From (3.19) in the proof of Proposition 
3.7 [8], 


y{t) < 4>{2f{t) + (1 + 8 C/ 37 ) 2-7 ) < 2 c(l + 2f{t) + (1 + /3) 2-7 [16(1 + c)] 2-7 y. 

1 

Estimating (1 + / 3 ) 2-7 by (jl.lll) . we have 

y{t) < 2c(l + 2f{t) + 25W-\l + /35^)[16(1 + c)]'^y < 2c[32(l + c)]5W(/(t) + 1 + p^y 

2(27+1) 2.7 ^27 

< 2 2-7 (1 + c)~ y + (I~ + fit)'^) < 3 • [2^(1 + c)]~ (1 + + /(t)"^), 


which proves (lA.ll) . 


□ 
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Lemma A.3. Let T G M. Suppose the continuous functions y{t),f{t) > 0 and h{t),g{t) > 0 on 
[T, oo) satisfy 

y'{t) < -h{t)y{t) + f{t) Vt > T, 

a'(t) 

/ hMdr = oo and lim = 0, (A.2) 

then 

limsup(g(t)y(t)) < lim sup ^ ) • 

i->-oo t^oo ^ n{t) / 

Proof. Same as Lemma A.3 of [12]. □ 

Lemma A.4. Let f{t) > 0 be a C^-function on (0,oo). Assume 

fd = limsup[/'(t)]“ < oo. 

t—>-oo 

Then there is T > 0 such that for any t 2 > ti> T, 

/(^i) < 7 (^ 2 ) + (^2 — ii)(/3 + !)• (A.3) 

Proof. There exists T > 0 such that for all t > T one has —f'{t) < /3 + 1. Let t 2 > ti> T. Then 

/(^i) = 7(^2) - [ f'{r)dT < f{t2) + [ (/3 + l)dr = f{t2) + (t2 - h){j 3 + 1 ), 

Jti Jti 

which proves ()A.3I1 . □ 
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